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A  new  structure  of  passive  direct  methanol  fuel  cell  (DMFC)  with  two  methanol  reservoirs  separated 
by  a  porous  medium  layer  is  designed  and  a  corresponding  mathematical  model  is  presented.  The  new 
designed  passive  DMFC  can  be  directly  fed  with  highly  concentrated  methanol  solution  or  neat  methanol. 
The  porosity  (epr)  of  the  porous  medium  layer  is  optimized  using  the  proposed  model.  Some  operation 
parameters  are  also  optimized  by  both  the  numerical  calculation  and  experimental  measurement.  The 
new  designed  DMFC  can  be  continuously  operated  for  about  4.5  times  longer  than  a  conventional  passive 
DMFC  with  the  optimum  parameters.  The  methanol  crossover  during  the  same  discharging  is  only  about 
50%  higher. 
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1.  Introduction 

Direct  methanol  fuel  cells  (DMFCs)  have  attracted  worldwide 
attention  due  to  their  high  energy  density  and  low-pollution  [1]. 
There  are  two  distinctions  for  fuel  delivery  system  for  DMFCs: 
active  and  passive.  In  the  active  DMFC,  both  methanol  solution  and 
air  were  supplied  to  the  catalyst  layer  through  flow  fields  using  a 
liquid  syringe  pump  and  compressed  gas  cylinder.  Compared  with 
the  active  DMFC  systems,  both  cathode  and  anode  flow  fields  were 
removed  in  the  passive  DMFC  systems  and  a  fuel  reservoir  was 
placed  at  anode,  methanol  was  delivered  from  fuel  reservoir  to 
anode  catalyst  layer  by  diffusion.  Therefore,  the  passive  DMFCs  are 
considered  as  the  most  promising  fuel  cells  for  the  portable  elec¬ 
tronic  devices  because  of  their  simple  structure  [2].  However,  there 
are  a  number  of  issues  have  to  be  solved  before  the  commercial¬ 
ization  of  DMFCs,  particularly  the  methanol  crossover  through  the 
proton  exchange  membrane  (PEM)  [3].  Plenty  of  works  have  been 
carried  out  to  reduce  methanol  crossover  [4-10]. 

DMFC  fed  with  dilute  methanol  solution  is  a  facile  way  to 
alleviate  the  problem  of  methanol  crossover  (i.e.,  1. 0-2.0  mol  L-1 
for  active  methanol  supply  systems  and  3.0-5.0molL-1  for  pas¬ 
sive  methanol  supply  systems)  [11-15].  However,  using  diluted 
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methanol  solution  as  the  fuel  for  passive  DMFCs  will  result  in 
a  short  operation  time  as  well  as  a  low  specific  energy  density 
because  of  the  constant  volume  of  the  fuel  reservoir.  Hence  it  is 
urgent  to  find  an  effective  fuel  feeding  strategy  which  is  suitable 
for  highly  concentrated  methanol  feeding  for  the  passive  DMFCs. 
Many  strategies  have  been  developed  to  increase  the  transfer  resis¬ 
tance  of  methanol  from  the  fuel  reservoir  to  anode  catalyst  layer 
(CL)  in  a  passive  DMFC  to  meet  the  objective  of  high  concentra¬ 
tion  methanol  feeding  [16-20].  Abdelkareem  et  al.  [16-18]  added 
a  dense  porous  carbon  plate  between  fuel  reservoir  and  anode 
CL  of  a  passive  DMFC.  The  data  showed  that  the  increased  mass- 
transfer  resistance  of  methanol  resulting  from  the  porous  carbon 
plate  makes  it  possible  to  feed  high  concentration  methanol  for  the 
passive  DMFC.  When  employing  a  16.0  mol  L-1  methanol  solution, 
a  maximum  output  power  density  of  24  mW  cm-2  can  be  yielded  at 
room  temperature.  In  contrast  to  the  conventional  passive  DMFCs, 
the  current  density  decreased  much  more  slowly  during  the  5.0  h 
discharging  time  under  0.1  V  cell  voltage.  Kim  et  al.  [19,20]  added 
hydrogels  into  the  fuel  reservoir  to  reduce  methanol  diffusion 
from  fuel  reservoir  to  anode  CL,  so  that  the  air-breathing  passive 
DMFC  can  be  operated  at  a  high  methanol  concentration  up  to 
8.0molL_1.  When  the  methanol  concentration  was  changed  from 
4.0  mol  L-1  to  8.0  mol  L-1 ,  the  maximum  power  density  of  the  pas¬ 
sive  DMFC  equipped  with  hydrogel  fuel  reservoir  was  increased 
from  1 6.5  mW  cm-2  to  21 .5  mW  cm-2,  while  the  maximum  power 
density  of  a  conventional  DMFC  decreased  from  14.0mWcm-2  to 
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Nomenclature 

rMEOH 

L0 

methanol  concentrations  in  DMR 

rMEOH 

Lac 

methanol  concentration  in  anode  catalyst  layer 

rMEOH 

Lad 

local  methanol  concentration  in  anode  diffusion 
layer 

Cc°2 

oxygen  concentration  in  cathode  catalyst  layer 

Cc°<f 

local  oxygen  concentration  in  cathode  catalyst  layer 

rMEOH 

Lh 

methanol  concentrations  in  HMR 

rMEOH 

Lpr 

local  methanol  concentration  in  porous  medium 
layer 

rMEOH 

Lref 

reference  methanol  concentration 

c% 

ref 

reference  oxygen  concentration 

rjMEOH 

ad 

diffusion  coefficient  of  methanol  in  anode  gas  diffu¬ 
sion  layer 

DS? 

diffusion  coefficient  of  oxygen  in  cathode  gas  diffu¬ 
sion  layer 

rjMEOH 

h2o 

diffusion  coefficient  of  methanol  in  water 

dmeoh 

diffusion  coefficient  of  methanol  in  PEM 

dme°h 

diffusion  coefficient  of  methanol  in  porous  medium 
layer 

E° 

cell 

ideal  electromotive  force  under  standard  condition 

F 

Faraday’s  constant 

j 

current  density  of  the  fuel  cell 

Ja 

current  density  in  anode  side 

Jet 

total  current  density  in  cathode  side 

/MEOH 

3rpf 

reference  exchange  current  density  of  anode 

& 

reference  exchange  current  density  of  cathode 

h 

thickness  of  reservoir  for  low  concentration 
methanol 

k 

thickness  of  reservoir  for  high  concentration 
methanol 

Nad 

methanol  flux  through  the  anode  diffusion  layer 

Ncd 

oxygen  flux  through  the  cathode  diffusion  layer 

Nm 

methanol  crossover  flux  through  PEM 

Npr 

methanol  flux  through  PML 

„h2° 

electro-osmotic  drag  coefficient  of  water 

„MEOH 

“d 

electro-osmotic  drag  coefficient  of  methanol 

R 

universal  gas  constant 

Rcon 

contact  resistance  of  cell 

T 

cell  temperature 

Greeks  letters 

anodic  transfer  coefficient 

ac 

cathodic  transfer  coefficient 

Ka 

order  of  anode  reaction 

y  c 

order  of  cathode  reaction 

$pr 

thickness  of  porous  layer 

thickness  of  membrane 

$ac 

thickness  of  anode  catalyst  layer 

<$cc 

thickness  of  cathode  catalyst  layer 

Sad 

thickness  of  anode  diffusion  layer 

Scd 

thickness  of  cathode  diffusion  layer 

G  m 

conductivity  of  PEM 

£ 

correction  factor  of  diffusion  coefficient 

£pr 

porosity  of  the  PML 

Peff 

fcpr 

effective  porosity  of  porous  medium  layer 

*7a 

anodic  overpotential 

cathodic  overpotential 

Air 


Cathode  DL 
Cathode  CL 
PEM 

Anode  CL 
Anode  DL 

DMR 

PML 

HMR 


Fig.  l.  Schematic  diagram  of  the  new  designed  passive  DMFC  with  two  fuel  reser¬ 
voirs.  (DMR:  diluted  methanol  reservoir;  PML:  porous  medium  layer;  HMR:  high 
concentration  methanol  reservoir). 


a  value  smaller  than  5.0mWcm-2.  A  semi-passive  DMFC  using 
porous  PTFE  plate  as  the  methanol  barrier  layer  (MBL)  was  manu¬ 
factured  by  Li  et  al.  [21,22].  A  maximum  output  power  density  of 
1 1 5.8  mW  cm-2  can  be  obtained  when  MBL  of  6.4  mm  and  Nation 
212  were  used  with  20  mol  L-1  methanol  supplied. 

Above,  all  the  previously  reported  designs  can  be  operated 
with  relative  high  concentration  methanol  solution.  However,  the 
employment  of  methanol-transfer-control  layer  will  bring  a  much 
larger  C02  releasing  resistance  than  conventional  passive  DMFC, 
because  C02  must  transported  through  the  same  components  as 
methanol  in  the  above-mentioned  designs.  In  order  to  realize  the 
passive  DMFC  supplied  with  high  methanol  concentration,  a  pas¬ 
sive  DMFC  with  two  methanol  reservoirs  separated  by  a  porous 
medium  layer  (PML)  is  designed  in  the  present  paper  as  shown  in 
Fig.  1.  The  methanol  crossover  can  be  well  controlled  in  this  DMFC, 
because  the  methanol  concentration  in  anode  catalyst  layer  can 
be  controlled  at  a  low  level  during  the  whole  discharging  process. 
Moreover,  C02  transport  in  this  DMFC  does  not  pass  through  the 
PML,  so  that  C02  releasing  resistance  in  this  new  designed  DMFC 
can  keep  the  same  as  in  the  conventional  DMFC.  The  continuous 
operation  performance  of  this  new  DMFC  is  investigated  by  both 
mathematical  model  and  experiments.  Compared  with  the  con¬ 
ventional  DMFCs,  the  new  designed  DMFC  can  get  much  better 
continuous  operation  performance.  The  structural  and  operation 
parameters  are  also  optimized  by  the  mathematical  model  in  this 
study. 

2.  Experimental 

PtRu/C  (60wt.%  Pt,  Pt:Ru  =  l:l)  powder  and  Pt/C  (60wt.%  Pt) 
powder  were  suspended  respectively  in  10%  PTFE  solution  in  an 
ultrasonic  bath  until  a  homogeneous  ink  was  formed.  After  spraying 
the  catalyst  ink  to  20%  PTFE  wet-proofed  two  2  cm  x  2  cm  carbon 
papers,  it  was  then  cured  at  340  °C  in  a  nitrogen-filled  vacuum 
oven  for  1  h  to  form  the  anodic  and  cathode  electrodes.  The  load¬ 
ings  of  noble  metal  and  Nation®  in  both  catalyst  layers  were  about 
4  mg  cm-2.  The  electrodes  were  hot  pressed  to  the  membrane  at 
130  °C  and  5  MPa  for  3  min  to  get  a  MEA. 

Polycarbonate  was  chosen  as  the  reservoir  material,  and  316L 
stainless  steel  palates  with  gold  layer  electroplated  on  the  surface 
were  used  as  current  collectors.  The  conventional  passive  DMFC 
was  fabricated  following  Ref.  [23],  the  dimensions  of  the  reservoir 
were  2  cm  x  2  cm  x  1  cm  (thick).  The  new  designed  DMFC  has  two 
fuel  reservoirs  with  the  dimensions  of  2  cm  x  2  cm  x  0.1 8  cm  (thick) 
and  2  cm  x  2  cm  x  0.72  cm  (thick)  as  shown  in  Fig.  1 ,  and  these  two 
reservoirs  were  separated  by  porous  medium  layer  (0.1  cm  thick 
PTFE  porous  membrane). 
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After  the  MEA  installed  in  the  cell,  methanol  solution  was 
directly  fed  into  the  fuel  reservoirs.  Both  the  polarization  curves 
and  discharging  curves  were  measured  at  25  °C. 

3.  Model  description 

3.1.  Assumptions 

The  schematic  diagram  of  the  new  designed  passive  DMFC  is 
shown  in  Fig.  1,  and  this  DMFC  is  named  as  H-DMFC  because  it  can 
be  fed  with  highly  concentrated  methanol  solution. 

The  assumptions  adopted  for  the  model  in  this  study  are  as 
follows: 

1)  The  operation  temperature  of  the  DMFC  is  assumed  to  be  uni¬ 
form  and  the  electrochemical  reaction  occurs  at  a  constant 
temperature  of  25  °C. 

2)  The  concentration  gradient  of  the  reactants  across  the  catalyst 
layers  is  ignored  because  the  thickness  of  the  catalyst  layer  is 
much  smaller  than  that  of  the  diffusion  layer. 

3)  The  methanol  which  permeates  from  anode  CL  to  cathode  CL  is 
assumed  to  be  completely  consumed  by  the  electrooxidation 
reaction  [24]  at  the  boundary  between  PEM  and  cathode  CL, 
where  consequently  the  methanol  concentration  drops  to  zero. 

4)  Methanol  concentrations  in  the  high  concentration  methanol 
reservoir  (HMR)  and  the  diluted  methanol  reservoir  (DMR)  are 
assumed  to  be  uniform,  respectively.  The  H-DMFC  discussed  in 
this  work  has  a  constant  HMR/DMR  volume  ratio  of  4:1. 

3.2.  Mass  transfer  equations 


In  the  cathode  diffusion  layer,  oxygen  flux  can  be  also  expressed 
by  Fields  law  because  of  no  external  convection  exists  in  the  cathode 
side  of  an  passive  DMFC,  so 

Wcd  =D°^  (5) 

where  Ncd  is  the  methanol  flux,  c°d2  is  the  oxygen  molar  concen¬ 
tration,  and  D°2  is  the  effective  diffusivity  of  oxygen  in  cathode 
diffusion  layer. 

3.3.  Conservation  equations 

The  methanol  flux  through  the  anode  diffusion  layer  is  con¬ 
sumed  through  two  routes:  methanol  electrooxidation  in  anode 
catalyst  layer  and  methanol  permeation  across  the  PEM.  The  elec¬ 
trochemical  reaction  in  anode  CL  of  DMFC  is 

CH3OH  +  H20  -*  C02  +  6H+  +  6e“  (6) 

According  to  this  reaction  formula,  1  mol  methanol  can  produce 
6  mol  protons  and  6  mol  electrons.  Thus 

Nad  =  TJ  +  Nm  (7) 

In  the  cathode  CL  of  DMFC,  two  reactions  are  proceeding  simul¬ 
taneously.  The  expected  reaction  is  the  oxygen  electroreduction, 
every  1.5  mol  02  reacts  with  6  mol  protons  and  6  mol  electrons 
through  the  following  equation: 


Methanol  permeation  flux  through  PML,  which  is  the  methanol- 
transfer-control  layer  in  H-DMFC,  can  be  described  by  Fields  law 


NPr  =  —D 


MEOH 

pr 


acgJEOH 

dx 


(1) 


where  Npr  is  the  methanol  flux,  c^}E0H  is  the  methanol  concen¬ 
tration,  x  is  the  coordinate,  and  D^rE0H  is  the  effective  diffusion 
coefficient  of  methanol  in  the  PML  and  can  be  estimated  by  Brugge- 
man’s  relation  [25]  as, 


nMEOH  _  c1.5nMEOH 
upr  ~  ^pr  ^H20 

where  spr  is  the  porosity  of  the  PML  and  D^EE)H  is  the  diffusivity  of 
methanol  in  liquid  water. 

Methanol  flux  across  the  anode  diffusion  layer  Nad  is  calculated 
by  the  following  equation 


1 .502  +  6H+  +  6e“  -*  3H20  (8) 

The  other  reaction  occurs  in  cathode  CL  is  the  consumption  of 
permeated  methanol  from  anode  side, 

f  CH3OH  +  H20  ->  C02  +  6H+  +  6e- 
\  1 .502  +  6H+  +  6e-  -*  3H20 

Therefore,  the  oxygen  flux  across  the  cathode  gas  diffusion  layer 
can  be  determined  by  the  conservation  equation: 

JvCd  =  T/+!jvm  (io) 


3.4.  Potential  equations 


Nad  = 


QrMEOH 
nMEOH  a  ad 

dx 


0) 


where  c^E0H  is  the  methanol  concentration,  and  L>^E0H  is  the  effec¬ 
tive  diffusivity  of  methanol  in  anode  diffusion  layer. 

The  methanol  crossover  through  the  PEM  is  caused  by  both 
methanol  diffusion  and  electro-osmotic  drag,  and  can  be  estimated 
as 


Nm  =  -D 


MEOH  ULm 


acME°H 


dx 


+  n 


meohJ 


(4) 


Butler-Volmer  expression  is  used  to  describe  the  relationship 
between  current  density  and  overpotential  in  both  anode  and  cath¬ 
ode  catalyst  layers 


HD 

exp(^,£)  <i2) 


where  cf^E0H  is  the  methanol  concentration,  D]^eoh  is  the  effec¬ 
tive  diffusion  coefficient  of  methanol  in  PEM,  Nm  is  the  methanol 
crossover  flux,  n^E0H  is  the  electro-osmotic  drag  J  is  operation  cur¬ 
rent  density  of  the  DMFC,  F  is  Faraday’s  constant.  First  term  on  the 
right-hand  side  of  Eq.  (4)  represents  the  methanol  diffusion  across 
the  PEM,  and  the  second  term  describes  the  methanol  permeation 
through  the  PEM  caused  by  the  electro-osmotic  drag. 


where  qa  and  qc  are  overpotentials  in  anode  and  cathode  side, 
respectively.  Ja  is  the  current  density  in  anode  side  and  Ja  =  J,  while 
Jet  is  defined  as  the  total  current  density  in  cathode  side.  The  per¬ 
meated  methanol  is  considered  to  be  consumed  by  electrooxidation 
reaction  in  the  cathode  catalyst  layer,  so 

Jct=J  +  6FNm 


(13) 
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Table  1 

Parameter  values. 


Parameter/symbol  (unit) 

Value 

Ref. 

Cell  temperature/T  (K) 

298.15 

Reference  methanol  concentration /c^E0H  (mol cm-3) 

1  x  10“3 

[26] 

Reference  oxygen  concentration/^  (mol cm-3) 

1  x  10“6  x  101325/(RT) 

[26] 

Diffusion  coefficient  of  methanol  in  anode  gas  diffusion  layer/D^E0H  (cm2  s-1 ) 

,d-°h 

[26] 

Diffusion  coefficient  of  methanol  in  water/D^E°H  (cm2  s-1 ) 

2.8  x  10_6exp(2436/353)  -  (2436/T)) 

[26] 

Correction  factor  of  diffusion  coefficient/^ 

0.55 

[27] 

Diffusion  coefficient  of  methanol  in  membrane/DjJE0H  (cm2  s-1 ) 

4.9  x  10-6exp(2436/333)- (2436/T) 

[26] 

Ideal  electromotive  force  under  standard  condition/E°dl  (V) 

1.213 

[28] 

Faraday’s  constant /F 

96,487 

Reference  exchange  current  density  of  anode/j^E0H  (A cm-2) 

0.011  x5ac 

[29] 

Thickness  of  cathode  catalyst  layer /<5ac  (cm) 

0.005 

Assumed 

Reference  exchange  current  density  of  cathode/j°2  (A cm-2) 

0.011  x<5cc 

[29] 

Thickness  of  cathode  catalyst  layer/5cc  (cm) 

0.005 

Assumed 

Thickness  of  reservoir  for  low  concentration  methanol//0  (cm) 

0.18 

Thickness  of  reservoir  for  high  concentration  methanol/^  (cm) 

0.72 

Electro-osmotic  drag  coefficient  of  methanol/nJf1E0H 

xMEOHnH20 

[30] 

Electro-osmotic  drag  coefficient  of  water/njb0 

2.9exp(1029/333)(1029/T) 

[27] 

Contact  resistance  of  cell/RCOn  (fi  cm2) 

0.8 

Assumed 

Conductivity  of  PEM/crm  (S  cm-1 ) 

0.073exp(1268((l/298)  —  (1/T))) 

[26] 

Order  of  anode  reaction/ya 

1 

Assumed 

Order  of  cathode  reaction /yc 

1 

Assumed 

Anodic  transfer  coefficient/aa 

0.5 

Assumed 

Cathodic  transfer  coefficient/ofe 

0.5 

Assumed 

Thickness  of  porous  layer /<5pr  (cm) 

0.1 

Thickness  of  rnernbrane/<5m  (cm) 

0.015 

Thickness  of  anode  diffusion  layered  (cm) 

0.003 

Thickness  of  cathode  diffusion  layer/<5cd  (cm) 

0.003 

3.5.  Numerical  analysis 


According  to  assumptions  3  and  4  in  Section  3.1 ,  Eqs.  ( 1  )-(4)  can 
be  rewritten  as  follows 


rMEOH  _  rMEOH 

Npr  =  D™  OHC_h - ^0 - 


°pr 


Nad=^M/°H- 


Nm=D^EOH^c - +n 


MEOH-Za 

d  p 


Ncd  =  ^|(ca0^b-c°?) 


(14) 


cmeoh  an(^  cmeoh  here  represent  the  methanol  concentrations  in 
HMR  and  DMR,  respectively,  and  c^b  is  the  oxygen  concentration 
in  the  ambient  air.  <$pr,  8ad  8m  and  <5cd  are  the  thicknesses  of  the 
corresponding  layers,  respectively. 

The  methanol  consumption  only  happens  in  the  anode  and  cath¬ 
ode  catalyst  layer.  Therefore,  when  the  fuel  cell  is  working  under 
a  constant  current  density,  values  of  cJ^E0H  and  c^E0H  during  the 
discharge  process  can  be  calculated  by  the  following  differential 
equations: 


0CMEOH 

dt 


N, 


pr 


9c0MEOH  Npr  -  Nad 

dt  In 


(15) 


The  initial  conditions  for  Eq.  (15):  c^-°  =  c1^  and  c^~°  =  cjjnt. 
Combining  with  Eq.  (14),  the  solutions  for  cJ^E0H  and  c^E0H  can 
be  easily  obtained. 

The  cell  voltage  is  determined  by  the  following  equation 


V  —  Ece  ii  —  TJa  —  rjc  —JRcon 


(16) 


where  Ecen  is  the  electromotive  force,  which  can  be  calculated  by 
Eceii  =  E°eu  +  (T  -  To)(dEce\\ldT),  and  E°eU  is  the  electromotive  force 
under  standard  conditions. 

4.  Results  and  discussion 

The  conventional  DMFC  investigated  in  this  work  has  a  fuel 
reservoir  with  a  thickness  of  1.0  cm,  while  the  thickness  of  DMR 
and  HMR  in  H-DMFC  is  0.1 8  cm  and  0.72  cm,  respectively.  The  cross- 
sectional  area  of  the  reservoirs  is  equivalent  with  the  active  area  of 
the  electrode. 

4.1.  Estimation  of  model  parameters 

Most  of  the  parameters  in  the  analytical  solution  could  be  found 
in  the  literatures  as  given  in  Table  1,  but  some  need  to  be  assumed 
to  fit  the  experimental  results.  Therefore,  the  comparison  of  the 
experimental  polarization  data  and  modeling  results  for  a  conven¬ 
tional  air-breathing  passive  DMFC  is  illustrated  in  Fig.  2a  to  estimate 
the  assumed  parameters.  Three  different  concentrations  are  fed  for 
the  DMFC.  The  calculated  results  agree  well  with  the  experimental 
data  except  at  the  region  of  activation  polarization. 

The  model  parameters  are  also  estimated  by  the  transient  cell 
voltage  test  under  constant  current  density  of  both  conventional 
DMFC  and  H-DMFC.  A  0.1  cm  thick  porous  medium  layer  (PML) 
is  used  for  the  H-DMFC,  and  the  porosity  of  this  PML  is  0.4.  The 
DMR  and  HMR  are  fed  with  2.0  mol  L-1  methanol  solution  and 
neat  methanol,  respectively.  The  conventional  DMFC  is  fed  with 
3.0  mol  L-1  methanol  solution.  The  operation  current  density  for 
both  the  H-DMFC  and  conventional  DMFC  is  60  mA cm-2.  It  is 
shown  in  Fig.  2b  that  the  experimental  and  calculated  results  have 
good  agreements  for  these  two  DMFCs.  However,  the  experimental 
continuous  operation  time  is  about  2.0  h  shorter  than  the  calcu¬ 
lated  one  for  H-DMFC.  This  difference  is  caused  by  the  methanol 
evaporation,  which  is  not  taken  into  account  in  the  mathematical 
model,  during  the  discharging  process.  One  can  also  conclude  that 
the  H-DMFC  has  a  much  better  discharging  performance  than  the 
conventional  DMFC  from  Fig.  2b. 
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Fig.  2.  Experimental  and  calculated  (a)  polarization  data  for  a  conventional  DMFC; 
(b)  transient  cell  voltage  data  for  a  conventional  DMFC  and  a  H-DMFC. 

The  rest  part  of  this  section  will  focus  on  the  transient  behav¬ 
ior  analysis  and  parametric  study  of  H-DMFC.  The  dependence  of 
the  H-DMFC’s  performance  on  porosity  of  PML  is  investigated  in 
details,  while  the  thickness  of  the  PML  is  maintained  at  0.1  cm. 
Meanwhile,  the  operation  parameters,  such  as  operation  current 
density,  methanol  concentration  are  also  discussed. 

4.2.  Effect  of  PML  porosity 

Porous  medium  layers  with  various  porosities  are  invited  in  the 
model  to  investigate  the  continuous  operation  performance  and 
to  further  optimize  the  porosity  of  porous  medium  layer.  Fig.  3 
gives  the  discharging  curves  of  H-DMFCs  with  PML  porosities  rang¬ 
ing  from  0.1  to  0.7  at  60  mA cm-2.  The  methanol  solutions  fed  in 
HMR  and  DMR  are  20.0  mol  L-1  and  3.0  mol  L-1,  respectively.  The 
discharging  time  increases  rapidly  with  the  increasing  of  the  PML 
porosity  (£pr),  and  gets  a  maximum  value  at  the  porosity  0.4.  When 
£pr  is  0.1,  the  H-DMFC  can  be  continuously  operated  for  only  about 
1.0  h,  which  is  even  much  shorter  than  the  discharging  time  of  the 
conventional  DMFC.  The  reason  is  that  the  methanol  supplement 
from  HMR  to  DMR  cannot  satisfy  the  methanol  consumption  by 
electrochemical  reaction  in  anode  CL.  When  £pr  comes  to  0.4,  the 
continuous  operation  time  increases  to  about  17h,  which  is  more 
than  3  times  longer  than  that  of  the  conventional  DMFC.  Further 
increasing  the  MPL  porosity  will  gradually  decrease  the  discharging 
time. 


Fig.  3.  Calculated  discharging  curves  of  the  H-DMFC  with  different  PML  porosities 
ranging  from  0.1  to  0.7  at  60  mA cm-2. 
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Fig.  4.  Methanol  crossover  fluxes  of  the  H-DMFC  during  the  discharging  process 
with  different  MPL  porosities. 

Fig.  4  compares  the  variation  of  the  methanol  crossover  flux 
with  time  at  different  MPL  porosities  during  the  discharging  pro¬ 
cess.  When  £pr  is  0.1,  methanol  crossover  flux  falls  rapidly  with 
time  because  of  the  low  methanol  permeation  velocity  through 
MPL.  Methanol  crossover  flux  decreases  with  time  but  increases 
with  MPL  porosity  when  £pr  is  smaller  than  0.3.  After  £pr  increases 
over  0.3,  the  methanol  crossover  flux  increases  during  the  first  hour 
because  of  the  quick  increasing  of  methanol  concentration  in  DMR. 
This  increment  is  caused  by  the  large  methanol  flux  through  MPL 
compared  with  methanol  crossover  flux.  The  methanol  crossover 
flux  gets  its  climax  at  about  1.0  h  and  then  decreases  with  time. 
When  £pr  is  0.7,  the  maximum  value  of  the  methanol  crossover  flux 
is  larger  than  4.0  x  1 0-7  mol  cm-2  s-1 ,  which  is  about  two  times  the 
quantity  of  the  initial  crossover  flux  value. 

To  optimize  the  porosity  of  PML,  average  power  density  dur¬ 
ing  the  discharging  process  is  calculated  and  shown  in  Fig.  5.  The 
case  of  £pr  =  0.1  can  be  eliminated  from  the  porosity  optimization 
because  of  its  ultra  short  discharging  time.  The  average  power  den¬ 
sity  keeps  growing  with  the  increasing  of  £pr  and  has  a  maximum 
value  when  £pr  =  0.5.  And  the  discharging  time  reaches  a  maximum 
value  at  £pr  =  0.4.  In  contrast  with  the  case  £pr  =  0.5,  the  discharging 
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Porosity 

Fig.  5.  Comparisons  of  average  power  density  and  continuous  operation  time  of  the 
H-DMFCs  using  different  PML  with  porosities  ranging  from  0.1  to  0.7. 


time  in  the  case  of  £pr  =  0.4  is  0.8  h  longer,  while  the  average  power 
density  is  only  0.1  mW  cm-2  smaller.  Therefore  0.4  is  considered  as 
the  optimum  value  for  the  PML  porosity. 

From  the  above  discussions,  one  can  conclude  from  the  calcu¬ 
lated  results  that  the  appropriate  using  of  PML  in  the  DMFC  can 
improve  the  discharging  time  notably  without  deteriorating  the 
output  power  density.  Taking  the  two  factors  into  account,  we 
choose  0.4  as  the  optimum  value  for  PML  porosity. 

4.3.  Effect  of  operation  parameters 

Operation  parameters  studied  in  this  work  contain  the  methanol 
concentration  fed  in  FIMR  and  DMR  and  the  operation  current  den¬ 
sity.  The  optimized  value  of  the  PML’s  porosity  is  used  in  this  section 
while  the  other  structural  parameters  keep  the  same  as  in  Section 
4.2. 

4.3 A.  Methanol  concentration  fed  in  HMR 

In  Fig.  6,  we  investigate  the  effects  of  methanol  concentration  in 
HMR  on  the  transient  performance  of  H-DMFC  by  both  experimen¬ 
tal  measurement  and  numerical  calculation.  The  DMR  is  fed  with 
3.0  mol  L-1  methanol  solution  and  the  operation  current  density  is 
60mA cm-2.  Fuel  fed  in  HMR  varies  from  10.0  mol L-1  methanol 
solution  to  neat  methanol.  The  discharging  curve  of  a  conventional 
DMFC  fed  with  3.0  mol  L-1  methanol  solution  is  also  shown  in  Fig.  6. 

The  experimental  results  agree  well  with  the  numerical  predi¬ 
cation  as  shown  in  Fig.  6.  This  result  can  further  verify  the  proposed 
mathematical  model.  The  discharging  time  increases  significantly 
with  methanol  concentration  in  HMR.  The  continuous  operation 
time  of  the  H-DMFC  has  about  5.0  h  increment  for  each  5.0  mol  L-1 
methanol  concentration  increasing  and  up  to  about  20.0  h  when 
neat  methanol  is  fed  in  HMR.  This  is  about  4.5  times  longer  than 
the  continuous  operation  time  of  the  conventional  DMFC  fed  with 
3.0  mol  L-1  methanol  solution. 

Meanwhile,  the  variations  of  methanol  crossover  fluxes  are 
shown  in  Fig.  7,  and  the  average  crossover  fluxes  are  given  in  the 
inset.  We  can  see  that  the  average  methanol  crossover  fluxes  in 
the  neat  methanol  case  is  about  200%  higher  than  the  flux  in  the 
10.0  mol  L-1  case,  while  it  is  only  about  50%  higher  than  that  of  a 
conventional  DMFC. 

Fig.  8a  is  the  bar  diagram  of  calculated  average  power  den¬ 
sity  and  continuous  operation  time  of  the  H-DMFC  with  different 
methanol  concentrations  fed  in  HMR.  Although  the  methanol 
crossover  increases  visibly  with  the  methanol  concentration  in 
HMR,  the  average  power  density  of  the  H-DMFC  does  not  deterio- 


Fig.  6.  Discharging  curves  of  a  conventional  DMFC  and  H-DMFC  with  different 
methanol  concentrations  fed  in  the  HMR  at  60  mA  m-2  (lines:  solid— 10  mol  L-1  (cal¬ 
culated):  dash— 1 5  mol  L_1  (calculated):  dot— 20  mol  L_1  (calculated);  dash  dot— neat 
methanol  (calculated);  symbols:  0—10  mol  L-1  (experimental):  A— 15  mol  L-1 
(experimental);  Q— 20molL_1  (experimental);  □— neat  methanol  (experimental)). 


Time(hour) 

Fig.  7.  Evolution  of  methanol  crossover  flux  for  a  conventional  DMFC  and  H-DMFC 
with  different  methanol  concentrations  fed  in  the  HMR  at  60  mAcrrr2. 


rate  with  the  increasing  of  methanol  concentration  because  of  the 
decreased  anodic  overpotential  as  shown  in  Fig.  8b.  It  can  be  seen 
from  Fig.  8a  that  the  average  power  density  even  slightly  increases 
with  the  methanol  concentration.  Therefore,  we  can  conclude  that 
neat  methanol  is  the  optimum  fuel  for  HMR,  because  of  the  longest 
discharging  time  as  well  as  the  highest  average  power  density  can 
be  obtained  in  this  case. 

4.3.2.  Methanol  concentration  fed  in  DMR 

The  calculated  discharging  curves  at  60 mAcm-2  with  differ¬ 
ent  methanol  concentrations  fed  in  DMR  and  neat  methanol  fed 
in  HMR  are  shown  in  Fig.  9a.  The  change  of  discharging  time 
with  the  methanol  concentration  in  DMR  is  neglectable  and  the 
five  discharging  curves  are  almost  identical  except  the  initial  3  h. 
The  variations  of  cell  voltage  during  the  first  3  h  are  enlarged  in 
Fig.  9b.  When  1.0  mol  L-1  methanol  solution  is  fed  in  DMR,  the  H- 
DMFC  has  the  lowest  cell  voltage  at  the  initial  stage  among  the 
five  cases  because  of  its  highest  anodic  overpotential  {pa)  as  shown 
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Fig.  8.  (a)  Bar  diagrams  of  average  power  density  and  continuous  operation  time;  (b) 
variation  of  anodic  overpotential  for  the  H-DMFC  fed  with  different  concentration 
methanol  solution  in  HMR. 


in  Fig.  9c.  After  0.5  h  discharging,  the  cell  voltage  for  1.0  mol  L-1 
case  increases  swiftly  and  becomes  the  highest  one,  because  the 
difference  of  among  the  five  cases  becomes  very  small  and  the 
cathode  overpotential  (rjc)  turns  to  be  the  determinant  factor  for 
the  difference  among  cell  voltages. 

Overall,  the  difference  among  transient  performances  is 
neglectable  when  the  methanol  concentration  in  DMR  is  changing 
from  1.0  to  5.0molL_1.  But  there  is  a  great  transient  perfor¬ 
mance  difference  in  terms  of  changing  the  concentration  in  HMR  as 
discussed  in  Section  4.3.1.  The  reason  is  that  changing  the  concen¬ 
tration  in  DMR  will  not  affect  the  methanol  concentrations  either 
in  DMR  (c0)  or  in  anode  catalyst  layer  (cac)  notably  during  the  dis¬ 
charging  process  as  shown  in  Fig.  10.  While  cac  is  the  determine 
factor  for  both  anodic  overpotential  and  methanol  crossover. 

The  methanol  crossover  fluxes  with  different  methanol  con¬ 
centrations  in  DMR  are  displayed  in  Fig.  11a  and  their  average 
values  are  plotted  in  Fig.  lib.  It  can  be  seen  from  Fig.  11a 
that  the  five  curves  cannot  be  differentiated  after  3.0  h  and  the 
average  value  of  methanol  crossover  flux  only  has  a  difference 
about  1 0-8  mol  cm-2  s-1  betweenl  .0  mol  L-1  and  5.0  mol  L-1  cases, 
which  is  caused  by  the  difference  of  methanol  crossover  fluxes  in 
the  first  3.0  h.  Meanwhile,  the  comparison  of  average  power  den¬ 
sities  during  the  discharge  process  is  also  shown  in  Fig.  lib.  When 
DMR  is  fed  with  2.0  mol  L-1  methanol  solution,  the  maximum  aver¬ 
age  power  density  is  obtained,  which  is  about  0.08  mW  cm-2  higher 
than  that  of  the  5.0  mol  L-1.  Therefore,  2.0  mol  L-1  is  chosen  as 
the  optimum  methanol  concentration  fed  in  DMR,  which  is  diluter 
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Fig.  9.  (a,  b)  Evolutions  of  cell  voltage  of  the  H-DMFC;  (c)  anodic  and  cathodic 
overpotential  with  different  methanol  concentrations  fed  in  DMR. 


than  the  commonly  used  ones  (3.0-5.0molL-1)  in  conventional 
DMFCs. 

4.3.3.  Operation  current  density 

The  methanol  crossover  flux  of  H-DMFC  depends  strongly  on  the 
operation  current  density,  because  both  the  methanol  consumption 
by  electrochemical  reaction  in  anode  CL  and  methanol  perme- 
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Fig.  10.  Evolutions  of  methanol  concentration  in  (a)  DMR;  (b)  anode  catalyst  layer 
with  different  methanol  concentrations  fed  in  DMR. 
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Fig.  11.  (a)  Variations  of  methanol  crossover  flux;  (b)  comparisons  of  the  average 
power  density  and  average  crossover  flux  during  the  discharge  process. 


ation  caused  by  electro-osmotic  drag  are  determined  by  current 
density.  So  the  changes  of  methanol  crossover  flux  of  H-DMFC  at 
different  current  densities  are  calculated  out  and  shown  in  Fig.  12. 
And  the  optimum  PML  porosity  and  the  methanol  concentrations 
in  HMR  and  DMR  are  employed  for  this  methanol  crossover  cal¬ 
culation.  Compared  with  the  40  mA  cm-2  case,  the  peak  value  of 


Fig.  12.  Evolution  of  methanol  crossover  flux  of  the  H-DMFC  at  different  current 
densities. 


Current  Density(mAcm‘2) 

Fig.  13.  (a)  Evolution  of  cell  voltage  (b)  average  power  density  and  total  electric 
energy  output  of  the  H-DMFC  at  different  operation  current  densities. 
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the  methanol  crossover  flux  for  100  mA  cm-2  case  is  about  50% 
smaller,  and  the  average  crossover  is  15%  smaller.  This  difference 
of  methanol  crossover  flux  is  caused  by  the  lower  methanol  con¬ 
centration  in  anode  CL  when  higher  current  density  is  used. 

The  corresponding  cell  voltage  variations  obtained  by  both 
experimental  measurement  and  numerical  calculation  are  shown 
in  Fig.  13a.  The  calculated  curves  show  very  good  agreements  with 
the  experimental  data  for  all  the  four  current  densities.  It  can  be 
found  that  the  higher  operation  current  density  is  used,  the  shorter 
discharging  time  and  lower  cell  voltage  can  be  obtained.  In  the  case 
of  40  mA  cm-2 ,  the  cell  voltage  of  H-DMFC  can  maintain  above  0.3  V 
for  more  than  20  h,  and  stop  at  about  25  h.  When  the  operation 
current  density  increases  to  1 00  mA  cm-2,  the  FI-DMFC  can  contin¬ 
uously  be  operated  for  more  than  10  h  with  the  cell  voltage  above 
0.1 5  V.  The  average  power  densities  and  electric  energy  outputs  per 
square  centimeter  at  different  current  densities  are  calculated  and 
displayed  in  Fig.  13b.  As  seen  from  the  figure,  the  smaller  current 
density  is  used,  the  higher  electric  energy  output  can  be  obtained. 
This  high  electric  energy  means  a  high  utilization  rate  of  fuel.  The 
reason  is  that  the  higher  current  density  is  used,  the  more  energy 
is  lost  as  waste  heat.  The  average  power  density  has  a  maximum 
value  of  16.4mWcm-2  at  80  mA  cm-2,  while  the  electric  energy 
output  keeps  decreasing  with  increasing  operation  current  den¬ 
sity.  The  decrease  of  electric  energy  output  is  primarily  caused  by 
the  difference  of  the  continuous  operation  time  for  these  cases. 

Taking  discharge  performance  into  account,  it  is  found  that 
the  FI-DMFC  with  optimized  structural  and  operation  parameters 
can  obtain  much  longer  discharging  time  than  the  conventional 
DMFC.  And  the  highest  power  density  can  be  obtained  at  about 
80  mA  cm-2. 

5.  Conclusions 

Porous  medium  layer  is  used  as  methanol-diffusion-control 
material  in  our  new  designed  passive  DMFC.  The  porous  medium 
layer  positioned  between  two  separated  fuel  reservoirs  makes  it 
possible  to  feed  high  concentration  methanol  solution  or  neat 
methanol  in  the  anode  side  of  the  DMFC  directly.  The  advantage 
of  the  new  design  is  verified  by  the  transient  performance  com¬ 
pared  with  a  conventional  passive  DMFC  at  the  same  conditions 
and  the  continuous  operation  time  is  5.5  times  the  length  of  the 
continuous  operation  time  of  a  conventional  DMFC. 

A  numerical  model  is  also  developed  and  a  parametric  study 
is  conducted  to  investigate  the  effects  of  structural  and  opera¬ 
tion  parameters  on  the  transient  performance  of  the  new  designed 
DMFC.  It  can  be  concluded  as  follows: 

1)  The  best  continuous  transient  performance  can  be  obtained 

when  the  porous  medium  layer  with  a  porosity  of  0.4  is  used 

in  H-DMFC. 

2)  The  transient  performance  of  the  H-DMFC  depends  on  the 

methanol  concentration  in  HMR  intensely.  While  the  effect  of 


the  methanol  concentration  in  DMR  on  the  performance  of  the 
H-DMFC  is  almost  neglectable.  Neat  methanol  and  is  consid¬ 
ered  to  be  the  optimum  fuel  for  HMR  by  the  verification  of  both 
numerical  calculation  and  experimental  measurement. 

3)  Highest  power  density  output  of  the  H-DMFC  can  be  obtained  at 
the  current  density  about  80  mA cm-2. 
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